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ABSTRACT 

We present the first measurement of the angular correlation function in a if -selected 
galaxy survey, from two fields covering almost 10 square degrees. The angular corre- 
lation function at small angles is well described by a 9~ Q S power law, as for optically 
selected samples. The clustering amplitude is reduced by a factor of ~ 4 between 
galaxies with K < 15 and a fainter magnitude slice with 15 < K < 16. This allows 
us to place constraints upon the redshift distribution of the galaxies and their spatial 
correlation function. We find no clear evidence for a change in clustering amplitude 
when galaxies are selected by their observed B — K colours. 

Key words: surveys-galaxies: clustering -dark matter - large-scale structure of Uni- 
verse - cosmology:observations - infrared:galaxies 



1 INTRODUCTION 

The measurement of angular correlations in optically se- 
lected samples of galaxies has been a useful probe of cosmo- 
logical models and theories of galaxy formation. The shape 
of the angular correlation function, w(9), on large scales 
measured in the APM Survey (Maddox et al 1990, 1996) 
showed that there is more structure in the galaxy distribu- 
tion on large scales than is predicted by the standard Cold 
Dark Matter model. The scaling of the correlation amplitude 
with magnitude in faint samples has given some indication 
of the nature of faint blue galaxies {e.g. Efstathiou et al. 
1991, Roche et al. 1993). 

It is now possible to survey large areas of sky in the 
near-infrared K band. Measurement of the angular correla- 
tion function allows us to place constraints upon the shape 
of the redshift distribution of Jf-selected galaxies and on 
their spatial two point correlation function. 



2 MEASURED CORRELATIONS 

We have measured the clustering of galaxies in two fields 
imaged in the K band with a NICMOS3 detector and in the 
B, V and I bands using a 2048 2 CCD camera. Details of 
the image reduction can be found in Gardner et al (1996 in 
preparation). The star/galaxy separation for objects with 
K < 15, I < 18, V < 19 and B < 20 was carried out 
using a combination of image morphology and colour se- 
lection (Gardner 1995), and each galaxy identification was 
confirmed by eye. In the range 15 < K < 16, the star/galaxy 
separation was performed using colour alone. The NEP field 
measures 1.97° x 2.58°, while the NGP field is rectangular 
with sides 6.91° x 0.66°. The NGP field contains a galaxy 



cluster, so we have also estimated the angular correlations in 
this field with the ~ 0.7 square degrees containing the clus- 
ter removed. This effectively produces two disjoint regions 
which we designate NGP-CL. 

The angular correlations were measured using three dif- 
ferent estimators, which we shall refer to as the direct Wd{9), 
Hamilton Wh{9) (Hamilton 1993) and ensemble w e {8) esti- 
mators. In order to compute the direct and Hamilton esti- 
mates, 20000 points were put down at random within the 
survey regions. The ensemble estimator is calculated by di- 
viding the fields into a grid of cells of side 0.02°. The form 
of the estimators is given by 



w d {6) 
w h {8) 
w e {6) 



2N R (DP) 
{N G - 1) {DR) 
A{DD){RR) _ 

{DR} 2 
{N.Nj) , 



{Ni){Nj) 



(1) 
(2) 
(3) 



where No is the number of galaxies in the field that satisfy 
the specified magnitude limits and Nr is the number of ran- 
dom points. The angular brackets indicate an average over 
all pairs of points in the sample for the case of Wd and Wh, 
and an average over pairs of cells for the case of w e , that 
are separated by angle 6 ± 89. The quantities DD, DR and 
RR are the number of distinct data-data, data-random and 
random-random pairs at angular separation 6. The number 
of galaxies in cell i is given by A^. 

The direct estimator suffers from linear biases arising 
from spurious correlations between the survey window and 
the data and is restricted by the accuracy with which the 
mean density of galaxies is known (Hamilton 1993, Landy 
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Table 1. Measured angular correlations: The amplitude of a power law fit w(0) = A9 8 is given along with the 2a range. The final 
column gives the value of the integral constraint. 



sample 


field 


number of 


A 


A (high) 


A (low) 


a 2 






galaxies 










K < 15 


NGP 


907 


0.043 


0.047 


0.039 


0.041 




NGP-CL 


653 


0.017 


0.021 


0.011 


0.016 




NEP 


675 


0.035 


0.039 


0.029 


0.042 


15 < K < 16 


NGP 


2641 


0.0096 


0.0108 


0.0088 


0.0093 




NGP-CL 


2144 


0.0081 


0.0093 


0.0069 


0.0081 




NEP 


2122 


0.0084 


0.0098 


0.0070 


0.0102 



& Szalay 1993). The Hamilton and ensemble estimators do 
not suffer from such a bias. The mean galaxy density is 
determined from the galaxy counts in each field, and so all 
the estimators will be biased low with respect to the true 
angular correlations by an amount 

a 2 » ^ J dnicHW(0i2), (4) 

where wt(9) is the true angular correlation function and 
O is the solid angle of the field (Groth & Peebles 1977). 
We estimate the size of this 'missing variance' or integral 
constraint by assuming that the true angular correlations 
are given by Wt(0) = A9~ ' 8 and performing a Monte-Carlo 
integration of equation (^) over each field. The amplitude A 
is then set by fitting the measured correlation function on 
small angular scales to (Roche et al 1993): 

w(0) =A6- 08 -a 2 . (5) 

We have also estimated the magnitude of the integral con- 
straint by calculating the variance in the number of galax- 
ies found in mock versions of our fields drawn from a 
much larger simulated APM catalogue (Gaztanaga & Baugh 
1996). The integral constraint found in this way is approxi- 
mately 25% smaller, which indicates the error resulting from 
the assumption that w(9) is a power law up to scales corre- 
sponding to the angular size of the field. 

The angular correlations for galaxies with K < 15 
are shown in Figure ^ for each field. Throughout this pa- 
per we plot the measured correlations plus the estimated 
value of the integral constraint, given in the final column 
of Table |l|. We have measured w(9) in logarithmic bins of 
width Alog# = 0.2. The dotted line shows the angular 
correlation between the random points in each field and 
the data which biases the direct estimator. The dashed 
line shows the power-law correlation function given by a 
least squares fit to the measured correlations on scales of 
0.005° - 0.06°, as listed in Table | for the Hamilton esti- 
mator. The angular correlation function bins are not inde- 
pendent, so we quote the 2a range in the fitted amplitude 
in columns 5 a nd 6 of Table Errorbars are computed us- 
ing 5w(6) = y/(l + w(6))/{DD) (Hewitt 1982). This may 
underestimate the errors, so we plot the 2a errors obtained 
from this formula. We have found that the 2a errors com- 
puted in this way are comparable to the errors obtained with 
a bootstrap resampling of the data (Ling, Frenk & Barrow 
1986). 

The correlations in the NGP region rise above the power 



law fitted to small scales in the angular range 0.06° — 0.50°. 
The NEP field shows a deficit of pairs on these scales, com- 
pared with the fitted power law. This suggests that this field 
could be slightly underdense, as a break in w(6) is not ex- 
pected at 9 < 1° for K < 15; in the optical, this corresponds 
to a magnitude limit of approximately B « 19, and w{9) is 
observed to be a power law until 9 w 2° (Maddox, Efstathiou 
& Sutherland 1996). The agreement between the values of 
w{9) measured with the different estimators is impressive 
in the NEP and NGP fields. There are some discrepancies, 
however, between the estimators in the NGP-CL field due 
to the hole caused by omitting the cluster, which results in 
less reliable measurements of w(0) in this field. From now on 
we shall plot only the angular correlation function obtained 
from the Hamilton estimator. 

The angular correlation function in the K < 15 sample 
is compared with that in the deeper 15 < K < 16 sam- 
ple in Figure ^, where we have plotted Wh(9), again with 
2a errors. The dashed lines show the power-law correlation 
function fitted to the measurements at small angles. There is 
a reduction in the amplitude of clustering by a factor of ~ 4 
in the 15 < K < 16 sample compared with the brighter sam- 
ple. The correlations in the deeper NEP and NGP samples 
are better described by a power-law on all scales, suggest- 
ing that the sampling fluctuations which affect the brighter 
catalogue are no longer significant for K < 16. 

Finally, we examine whether the clustering strength de- 
pends on galaxy colour. We find no consistent difference in 
the amplitude of w{9) in each field when we split the K < 16 
catalogue into 'red' and 'blue' galaxies at the median B — K 
colour. Figure |^ compares the angular correlation function 
in the blue and red subsamples with the correlations for 
the whole K < 16 sample. In the NGP field, blue galax- 
ies do appear to be more strongly clustered than both the 
red galaxies and the full sample. However, the angular cor- 
relations in this field are dominated by the presence of the 
cluster. In the region that contains the cluster, there are 
~ 2 times more galaxies with blue colours compared with 
the number that have red colours. In the NEP field, the red 
galaxies appear to be more strongly clustered than the blue 
galaxies, in particular around log 9 ~ —0.6. 



3 INTERPRETATION 

The angular correlation function is related to the spatial 
correlation function £(r) and the redshift distribution of the 
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Figure 1. The measured angular correlation functions with the 
estimated integral constraint added for the NEP, NGP and NGP- 
CL regions for galaxies brighter than K = 15. The open squares 
show w(8) calculated with the direct estimator, the filled squares 
with the ensemble estimator and the stars with Hamilton's es- 
timator. The error bars on the direct estimator values are 2<r 
errors calculated as described in the text. The dashed line shows 
the power-law fitted to the measured w(6) on small scales. The 
dotted line shows the random-galaxy correlations. The arrows in- 
dicate the smallest linear dimension of each field. 



Table 2. The correlation length needed to reproduce the am- 
plitude of clustering in the NEP field, for a range of theoretical 
rcdshift distributions and for two different cases of the evolution 
of clustering with rcdshift 



rcdshift 


-ft" < 15 


15 < 


K < 16 


distribution 


ro 


e 


ro 


e 


optical 


6.8 


0.0 


4.5 


0.0 


fit 


6.3 


-1.2 


4.0 


-1.2 


passive 


7.5 


0.0 


6.1 


0.0 


evolution 


6.9 


-1.2 


5.2 


-1.2 


no evolution 


6.3 


0.0 


4.3 


0.0 




5.8 


-1.2 


3.8 


-1.2 



Figure 2. Angular correlation functions for galaxies with K < 15 
(open squares) and 15 < K < 16 (filled squares), with 2a error 
bars calculated as described in the text. The dashed lines show the 
power law w(0) fitted at small angles. The angular correlations 
are computed using the Hamilton estimator. 



galaxies by Limber's equation (Peebles 1980, §56). On small 
scales, both w(9) and £(r) are measured to have power law 
forms (e.g. Davis & Peebles 1983, Maddox et al 1990, 1996, 
Baugh 1996) and Limber's equation can then be written in 
the form (Efstathiou et al. 1991): 

r(( 7 -l)/2) A ^ 



w(9) — \/tv 
where 



r( 7 /2) e-r 



-1 ' J 



Jo 



dz/ 



(djV/dz)dz 



(6) 



(7) 



and 

g{z) = (dz/dx)x 1 -''F(x)(l + z)- (3+e - 7) , (8) 

where F(x) depends upon the cosmology and is equal to 
unity for Q — 1. Here we have parameterised the redshift 
evolution of the spatial correlation function as 



"(3+€) 



(9) 



If clustering is fixed in proper coordinates then e = 0. If 
the pattern of clustering is fixed instead in comoving coor- 
dinates, then e « —1.2 for 7 w 1.8 (Efstathiou et al 1991). 
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Figure 3. Angular correlation functions for galaxies selected by 
B — K colour from the K < 16 sample. The stars show w(6) 
for all galaxies. The filled circles show w(9) for galaxies redder 
than the median B — K colour, while open circles show w(B) for 
galaxies bluer than the median B — K. The error bars are 2<r 
errors calculated as described in the text. 



The redshift distribution of i^-selected samples is less 
well determined than that of optical samples. In order to 
model the redshift distribution, we have combined redshifts 
from the Hawaii if-band Survey (Songaila et al 1994) and 
from Glazebrook et al. (1995), weighting each galaxy by the 
field area, sampling rate and redshift completeness fraction 
as described in Baugh & Efstathiou (1993). We have used 
the total K magnitudes listed in Table 3 of Songaila et al 
and the 20h~ kpc metric aperture magnitudes from Table 
Al of Glazebrook et al. Figure ^ shows the observed red- 
shift distribution for several cuts in K magnitude. The solid 
curve shows the prediction of the parametric form adopted 
by Baugh & Efstathiou (1993) as a fit to the redshift distri- 
bution of optically selected galaxies 

(AN/Az)Az = ^^"V exp(-(^/ 2c ) 3/2 )dz, (10) 

where the median redshift is z m = 1.412z c . For the com- 
bined 13 < K < 15 dataset, we find a median redshift 
of z m ~ 0.155, whilst for the deeper slice 13 < K < 16 
we obtain z m ~ 0.175. Also plotted are the predicted red- 
shift distributions for 'passive' (dotted line) and 'no evo- 
lution' (dashed line) models constructed with the Bruzual 



-a 4 




13 < K < 15 



■f t 1 1 f- 



13 < K < 16 




15 < K < 16 




Figure 4. The histograms show the redshift distribution com- 
bined from the surveys of Songaila et al. (1994) and Glazebrook 
et al. (1995). The solid curve shows the parametric fit described 
in the text. The dashed curve shows the prediction of a no evo- 
lution model and the dotted curve shows a model with passive 
evolution. The circles show the predictions of the semi-analytic 
model for galaxy formation of Cole et al (1994). 



and Chariot (1993, 1996 in preparation) spectral energy dis- 
tributions described in Gardner et al (1996) and Gardner 
(1996). The circles show the predictions of the ab initio semi- 
analytic model for galaxy formation described by Cole et al 
(1994). The theoretical curves give reasonable fits to the 
combined dataset, except perhaps at higher redshifts, where 
in any case the observed redshift distribution is less well 
determined. 

A further constraint on the form of the redshift distri- 
bution is the scaling of the amplitude of w(8) with apparent 
magnitude between the K < 15 and 15 < K < 16 samples. 
In Table ^ we find the values of the correlation length ro 
that reproduce the observed amplitude of the angular cor- 
relation function in the NEP field, for two specific cases of 
the evolution of clustering, e — and e = —1.2. Local opti- 
cally selected samples give ro ~ 5/i~ 1 Mpc, though the real 
space spatial correlation function £(r) has a shoulder fea- 
ture between 4 < r < 25/t -1 Mpc (Guzzo et al 1991, Baugh 
1996, Ratcliffe et al. 1996), which complicates fitting a sim- 
ple power law to £(r). 
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4 SUMMARY 

We have presented angular correlation functions, w(9), for 
galaxies selected in the A"— band over a large area of sky. We 
find that w(9) is well fitted by a #~ 0,8 power law on small 
angular scales, as observed for optically selected galaxies. 
There is a factor of « 4 reduction in the amplitude of clus- 
tering for a faint magnitude slice, 15 < K < 16, compared 
with galaxies with K < 15. This behaviour can be repro- 
duced assuming the redshift distribution predicted by a sim- 
ple no evolution model, if the correlation length of galaxies 
with a median redshift of z ~ 0.23 is around 4.3/i -1 Mpc 
whilst for galaxies at a lower median redshift of z ~ 0.13, 
ro ~ 6.3/i _1 Mpc, and clustering is fixed in proper coordi- 
nates. We find no consistent change in the strength of clus- 
tering in the different fields when galaxies are selected by 
their observed B — K colours, but this is does not rule out 
a dependence of clustering amplitude on colour in samples 
of higher median redshift (Landy el al. 1996). 

The measured clustering contributes to the uncertainty 
in the number counts (Gardner et al 1996) because different 
patches of sky contain different numbers of galaxies depend- 
ing upon whether they are in high or low density regions. 
The integral constraint a 2 gives an indication of the frac- 
tional variance in the number counts (see Peebles 1980 §36). 
Taking into account the fact that we have two widely sep- 
arated fields, the rms variance at K < 15 is ay ms = 0.14 
whilst for 15 < K < 16 &rms — 0.07. These values are up- 
per limits, though tests with simulated catalogues that have 
the same clustering show that these numbers would be re- 
duced by ~ 0.02. The contribution to the variance in the 
counts from galaxy clustering is several times larger than 
that expected from Poisson statistics alone. 
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